liiiiilliil^ 

iililiiiiiBiiili^^ 



liiiii^liiiiiilil 



111 



TH| UNSTEADY LIFT 
A WIN^'OF FINITE ASPECT RATIO 



By.I^OBERT T. JONE$ 




1940 



For sale by the SapeHntond^iit of Boenments, "Wtmhlngium, B. C. 

Sabscription priee* $3. per year 



Price 10 cento 



iiiiii 



■■11! 



iiiiiiiiiiil 



iiiPiiiiiiii ilill^^ 

AERONAUTIC SYMBOLS 





Symbol 


, — : — 

Metric 


SiiiiiiiiiiiiB 


Vnit 


Abbrevia- 
tion 


Unit 


. Abbrevia- 
tion 




1. -■ , 1 

F 




m 

Kg',.-- . 


foot (or mile) 

second (or hour)—..--. 

weight of 1 pound- 

a, }run>r','^",- ■ K 


ft. (or mi.) 
sec. (or hr.) 
lb. 




weigM of 1 kilogram 




\ P 


horsepower (metrlo) « 4 C ^- 


• m,j).s. 


horsepower I ^ ^ « - ^« 


^p.h. 

^^^^^ ^i^^ - 


\ipet0rs per secon^^ . a- ^ 


* feet^per seoo^C-.^^^^-^-^^ 



2. GENERAL SYMBOLS \ 

I-, Kinematio vjaoosity 




^andard-, aceeMatiM of fgj:ay%=«&.8066§ 

Masf'=-^, v" , ^ ' Specific weight of "ataodaJrd" ajr, i.«265/l£g/m» or 

\Moip.ettt of ia«t€ia=«mF. (Indicate axis of 0.07651 Ib./ou. ft. .":7; \ . . - ; 

■ radtoa of gyration & by proper subscript.) . ^ „ - , . , , . 

Coefficient of tiacosity • , 

, 3. AERODYNAMIC SYMBOLS 



Area y 

8pm . , 



6, 

X, Mft, absolute coefficitot Oi'*^ 

0, \ l>rag? .absolute coefficient 0x>=^^ 

; ^^pfiJ^ <i^^g^ aib^olate coefficient Od^^ 



Do 



Z>i, Induced drag, absolute coefficient Obc^^ 
jp^, ara^ite drag, absolute coefficient O^p^^ 

12, Resultant force ^ 



i^, . An0e of setting of w«]^ Cr^f^^ to thrust ;| 

line) / ^y^' i , . ] 

i^, ] Angle of stabili^0r setting C^^^ t^pr thrust ji 
{ ^^.^ lino) > ^ ' ^^^'^ ^ \ ^^'^ . . 
Resultant moment ^ ' ^ " ^ 
a, . Resultant angular Telocity ^ 

p2|/ Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100, 
in.p,h, normal. p:pssure at 15^ C, the cor- l 
l^^^onding mittib^is 204,000; or for a ^de| J 
' . ol^iO cm chord, 40 m,p,s., the correspO|kdin| I 
Jibber is 274,000) ^ > ' ^ \ - 

tMrter-of-pressure coefficient (ta1a<^ of d^^pi^^^^^^ 

of c.p. from leading edge to choi'd length) : 
Angle of attapk 
Angle of downwash 

Angle of attack, infinite aspect ratia . , 

Angle of attack, induced ' ^ .l . a ^ 

Angle of attack, absolute (itteaau%d ttmot zhm^ 
lift position) ' _ 

^^^^ . Flight-path angle 



a, 
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THE UNSTEADY LIFT OF A WING OF FINITE ASPECT RATIO 

By Robert T. Jones 



SUMMARY 

Unsfeady4ift functions for wings of finite aspect ratio 
have been calculated by correcting the aerodynamic inertia 
and the angle of attack of the infinite wing. The cahula- 
tions are based on the operational method. 

The starting lift of the finite wing is found to be only 
slightly less than that of the infinite wing; whereas the 
final lift may be considerably less. The theory indicates 
that the initial distribution of lift is similar to the final 
distribution. 

Curves showing the variation of lift after a sudden unit 
change in angle of attack^ during penetration of a sharp- 
edge gust, and during a continuous oscillation are given. 
Operational equivalents of these functions have been devised 
to facilitate the calculation of lift under various conditions 
of motion. As an application of these formulas, the 
vertical acceleration of a haded wing caused by penetrating 
a gust has been calculated, 

INTRODUCTION 

The two-dimensional potential theory of airfoils in 
nonuniform motion was given by Wagner (reference 1) 
and has been extended to problems involving the motion 
of hinged or flexible airfoils by Theodorsen (reference 2) 
and Ktissner (reference 3). 

In the case of steady motion, a correction is known to 
be necessary before the results of the two-dimensional 
theory can be applied to wings of finite aspect ratio. 
A theory for the unsteady lift of finite wings was devel- 
oped in reference 4. This theory has since been some- 
what improved mathematically by making use of 
operational methods in the solution of the integral 
equations. (See reference 5.) The present report 
combines this previous work md extends the fJieory to 
show the effects of gusts. 

THE INDICIAL LIFT 
INFLtiENCE OF THE WAKE 

Owing to the presence of circulation, the lifting wing 
leaves in its path a surface of discontinuity, the local 
vortex strength of which is determined by the rate of 
change of circulation taken both across the span and 
along the flight path. (See fig. 1.) The distribution 
of vorticity in the wake is determiued by the assump- 
tion that the flow field at each iostant conforms to the 



Kutta condition. An essential feature of the problem 
is the determination of the influence of this wake on 
the flow at the wing. 

It is important to note that the wake is supposed 
to remain plane and undistorted. As a consequence of 
this assumption, the effects of different wakes are 
additive, permitting the various flows to be built up by 
superposition. Thus, if the solution for the growth of 
the increment of lift following a sudden change of nor- 
mal velocity— or, what amounts to the same thing 
under the assumptions iavolved, a sudden change in 
the angle of attack — :is known, this solution may be 
used as the element in an integral that gives the lift in 
any variable motion. With this point in mind, atten- 
tion will at first be directed to the special case in which 
the wmg starts suddenly from rest at # =0 with the 




3farf 



FiGimia 1.— Flow <aused by wing startup from 



normal velocity w and the flight velocity Uq, the 
velocities remaiarug constant thereafter. 

LIFT NEAR THE START 

The starting lift of any wing may be expressed by a 
simple theorem based directly on the Kutta condition. 
As a consequence of this condition, the portion of the 
wake adjacent to the trailing edge must move as an 
impermeable extension of the wing surface. Thus, the 
first element of wake formed must move with the same 
normal velocity as the wing. The flow produced at 
the first instant is what ncdght be caused by thfe wing 
in process of growing wider at the rate Uq while moving 
downward with the velocity w. The starting lift may 
then be thought of as the reaction to uniform motion 
of the wing as a body with increasing mass: 

L=w^ (1) 



1 
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where m' is the mass representiog the aerodynamic 
inertia of the wing in normal motion. 

. In order to apply equation (1) to the finite wing, 
the inertia factor for such a wing must be known as a 
function of the width. Solutions for elliptic plates are 
given by the classical hydrodynamic theory, and these 
solutions can be used to represent approximately the 
initial rates of increase of inertia of wings of oval or 
elUptic plan form. 

The distribution of potential over each ohordwise 
section of an elliptic plate in normal motion has the 
game form as the corresponding two-dimensional 
potential. Thus 

(2) 



where E is the elliptic integral giving the ratio of the 
semiperimeter to the span. At the normal velocity 
w=E, the potential distribution over any chord is 
represented by a circle having the chord as diameter. 
(See fig. 2.) If the edge of the plate distorts into a 
slightly wider ellipse, the change in potential arising at 
any point will be measured by the difference between 
the original and the slightly expanded circles. (The 
change in the factor E during widening may be neg- 
lected for ordinary aspect ratios.) The pressure dif- 
ference across the plate with chai]^ing potential is given 
by the formula 



l/o 



Uo 




FiGTT&E 2.— The wake and the distribution of potential over the chord shortly after 

the start. 



For t(?=£? 



50 —X 

~= cot B 

and, from the geometry of the circle^ 

50 



(4) 



50 
5^ 



50 rfAc 



= CWr- = o cot (6) 



The pressiire across the plate with the normal velocity 
'w==^E and the flight velocity Z7o is, therefore, 



Vt^'^pUi^ cot cot 0 



(6) 



Integration of this pressure over any section gives the 
lift coefficient, for angle of attack a of the plate, 



(7) 



with each local center of pressure at the quarter-chord 
point. 

The start of the plane elliptic wing being equivalent 
to a uniform lengthening of each chord, the true elliptic 
outline is not preserved. Such a change, however, may 
be shown to conform very nearly to a change into 
another, slightly larger, ellipse at aU points except those 
very near the tips. Furthermore, if the wing is assumed 
to distort in any of a number of ways into a slightly 
different elliptical plan form, the change of aerody- 
namic inertia will be found to be but little affected by 
the change in shape and to depend primarily on the 
over-all change in size. Each such distortion can be 
thought of as representing a certain distribution of the 
starting velocity TJ around the edge of the wing. 
Equation (5) is ^act for all distortions of this class. 
Inasmuch as they may be made to fall on either side of 
the distortion represented by Z7= constant (represent- 
ing the start of the rigid wing), Uie equation is also 
conadered applicable to this case. 

THE DOWNWASH CORRECTION 

A reasonably accurate curve of the growth of Hft 
might now be drawn by connecting the starting value 




Basio flow 




Trarvsformed flow 

FiGUEE 3— Element of circulatory flow. 



(equation (7)) asymptotically to the known steady 
value given by the Prandtl theory. Calculations have 
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shown, however, that, after the wing has progressed a 
distance of the order of one semispan, the effect of 
finite width of the wake can be treated simplj, as a 
modification of the angle of attack of the entire wing, 
as in the steady-lift theory, A eloser approach to the 
true form of the curve may be obtained by proceeding 
on this basis. 

Before the three-dimensional problem is considered, 
it will be helpful to review certain aspects of the two- 
dimensional theory (reference 1). In order to make 
the analysis nondimensional, all velocities are expressed 
in terms of the flight velocity Uq and all lengths, in 
terms of the half chord. 

Figure 3 shows the elementary two-dimensional flow 
used as a starting point by Wagner (reference 1). 
This flow is caused by two vortices, representing, re- 
spectively, an element of circulation around the wing 
and the vortex left in the wake when this circulation 
originated, The streamlines of this flow are eccentric 
circles. One such circle (of unit radius) is chosen to 
represent the wing section and the axes are so placed 
that this circle has its center at the origin. The geom- 
etry of the resulting pattern is such that, when the 
wake vortex is at 0, the wing vortex will be at 
This spacing preserves the unit circle as a streamline 
of the flow. 

Transformation of the pattern by the formula 

2f=^-hi (8) 

flattens the unit circle into a thin-line wing section and 
distorts the originally circular streamlines into oval 
Joukowski figures. The transformed pattern thus rep- 
resents the circulatory flow around a flat wing section 
with an associated counter vortex in the wake. In the 
transformation, the centroid of wing vorticity remains 
at the position of the original bound vortex while the 
wake vortex is shifted forward somewhat as shown 
(fig. 3). 

Each elementary flow of the type shown contributes 
a certain velocity around the trailing edge of the airfoil. 
The flow due to an instantaneous change of angle of 
attack of the airfoil may be superposed on these flows 
and wiU contribute a trailing-edge velocity of opposite 
sense. On this basis, the problem of circulation with 
varying angle of attack may be solved by an inverse 
procedure. Assume some convenient distribution of 
wake vorticity and calculate (by integration) the trail- 
ing-edge velocity at each point along the flight path 
corresponding to the prescribed wake. The particular 
variation of angle of attack necessary to cancel this 
trailing-edge flow at each instant (Kutta condition) can 
then be determined. If a number of such curves are 
found, they may be added in various ratios so as to 
approximate some prescribed variation of angle of 
attack; the corresponding circulation curves are added 
in like ratios. 



In essentially the mancner described, Wagner (refer- 
ence 1) calculated the two-dimensional flow around a 
wing section foUowing a sudden unit change in angle of 
attack. The integrated pressures over the airfoil give* 
a lift coeiB&cient that as3iTnptotically approaches the 
known steady value 27r; whereas' the starting Hft 
coefficient is found to be exactly one-half this value. 
The center of pressure remains at the quarter-chord 
point throughout the motion. 

In the case of the finite wing, an element of the wake 
win be as depicted in figure 3 but will, in addition, 
contain vortices completing each circuit to the wing 
through the tips. The length of the tip vortices may 
be approximated by assuming that they extend to the 
chordwise centroid of the wing circulation. After some 



— I 




o w 



Figure 4.— Position of tlie centroid of discontinuity in the wing for different positions 
of the wake vortex. 

calculation, the equivalent length of the tip vortices 
in terms of the distance traveled s reduces to 

x==^^s{s+2) (9) 

Figure 4 illustrates the rapid travel of the centroid of 
discontinuity within the wing subsequent to its initial 
position at the trailing edge,^ It is seen that, after a 
travel of several chord lengths, the centroid may be 
taken at the middle of the wing section. This assump- 
tion will later be used. 

Figure 5 shows how an elementary loop vortex in the 
wake of a finite wing can be formed by cancelation from 
an element of the wake of an infinite wing. The 
downwash induced by segments CD and FH accounts 
for the aspect-ratio effect. Since a uniform distribu- 
tion of the downwash flow is assumed, the calculations 
will be restricted to the center of the wing. By the 
application of Biot-Savart's rule, the downwash velocity 
due to elements CD and FH is found to be 

^[(f ^i)vTO^"fl ^^^^ 

This expression for downwash may be integrated in 

1 At » ~0, the tip vortices are lengthening at an infinite rate and, although the vortex 
strength is zero at the beginning of the motion, the limiting calculation shows that the 
induced downwa^ flow has a certain rate of acceleration at this instant. As a result, 
the starting lift of the finite wing is diminished, in accordance with the result of the 
previous calculation. 



4 



iMBPO^lT NO. 681 — NATIONAL ADVISORY COiaMITTEE FOR AERONAUTICS 



B 

D 



< — ™~ . — ^ ^ 

FIG^RE 5.— Superposition of vortices to obtain finite loop. 

the case of elliptic spanwise loading. Let y represent 
the circulation aroiind any chordwise section; then 



7=r sin 6 



(11) 



where y~2 6j and F is the value of 7 at the center 
section. 

2/''+a:='=(|J[n-(5^J](l-We) (12) 
1 

where F?= 



1 + 



(6/2) 



Then the induced velocity due to a series of finite loops 
of the form CEF (fig. 5) is given by 



de 





6/2+ X """^ ^ IVF 



de 



-Psin^^ 



where X (A*) and jE (k) are the complete elliptic integrals. 
(See Peirce's table.) 

Subtracting the two-dimensional vortex E gives the 
effect of a series of segments of the form of DC and FH, 
distributed along the span according to the elliptic 
loading. 

+lE(k)^lJ (14) 

Figure 6 shows the variation of downwash velocity 
with increasing length of the wake as determined by 
this formula. Some additional rough calculations have 
shown that the downwash becomes practically uniform 
over the entire iving before the wake has attained a 
length of one semispan. 

Figure 7 shows downwash curves derived from equa- 
tion (14) for elliptic wings of aspect ratios A of 3 and 6. 
In this derivation, the unit of length was taken as half 
the central chord of the wing. Thus, the wings have 
the same chords (^^-^=2) but are of different spans. In 
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Figure 6.— 'Growth of downwash with increasing length of the wake. Pg— 1.0; 
elliptical span load. 

order to define the later portions of the curve, the wake 
was assumed to start with length equal to the mean 
chord hi A in each case. This assumption, though 
somewhat arbitrary, makes allowance for the curvature 
of the trailing edges of the wings. 




6 8 
s , tiglf ctiords 

Figure 7.— Downwash fujictioBS, «;j,(s). 

The iuduced downwash Wi with any variation of 
circulation r(s) along the flight path may be deter- 
mined from the curves given in figure 7 by superposition ; 
thus 

««?<«(«) ='W?r(s)r„CO) + J^ Wr(5~5o)r/(so)(^So (15) 

The growth of circulation following a sudden start of 
the motion will be determined from the two-dimen- 
sional theory by usmg the effective normal velocity 



(16) 



Let Fow be the rise of circulation followmg a sudden 
start with, unit normal velocity as given by the two- 
dimensional theory. Then, for the finite wing, 



r«,(s)=ro^,(s)~-ro^(s)i^,-(0)-- ro^(s-so)t^/(5o)c^5o (17) 

The determination of the effective normal velocity 
and the circulation for the finite wing thus depends on 
the simultaneous solution of integral equations (15) 
and (17). This solution may be conveniently obtained 
by operational methods. 
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OPERATIONAL SOLUTION OF INTEGRAL EQUATIONS 

Let D represent the operator djds and let 1=1 (s) 
represent the unit jump function, that is, a function of 
s having the yalue zero at s<^OVand haying the value 
1 at s>0. A f unction of s may be represented by a 
combination of operations on the unit jump function 



$(5)-i(Z))l 



(18) 



The combination of operations necessary to 

.reproduce the function ^{s) is called the operational 
equivalent of the function #(s). 

Kules for finding such equivalents are discussed in 
reference 6. The most general rule for proceeding 
either from # to or vice versa, is: 



The rule needed in the foUownig development is 
the Heaviside expansion theorem: 



where / and F are algebraic polynomials and the Vs 
are the roots of jF(X) =0. 

The operational treatment of integral equations is 
based on the proposition that an integral of the form 
of (15) may be regarded as the Knear superppsition of 
the effects of a succession of small jump functions. 
The operational form of (15) is 

'^UD)==Wr{D)fM . (20) 

and that of (17) 

f^{D):^fo^iD)[l~wUD)] (21) 
Solving algebraically for Wiy, (D): 



(22) 



The induced velocity Wi^{s) gives the variation of 
the effective angle of attack of the finite wing when the 
geometric angle of attack is held constant. The lift at 
later stages of the motion is then found by combining 
the effective angle-of-attack variation 



(23) 



with the two-dimensional indicial-lift function given 
by Wagner. Let Clo^(s) = (7lq^(s) be the lift in two- 
dimensional flow following a sudden unit jump of angle 
(the curve given by Wagner is for a=l/2w)] then, for 
the finite wing: 

^?z,^(«)=Gco«(^) - Ci,o,(s)Wi„(0) ~JjGL^,is-'So)iVi^'iso)dso 

(24) 



or, in operational form: 



(25) 



Substitution of the expr^sion ior WuoiP) froin equation 
(22) gives the operational form of the lift function for 
the finite wing in terms of the known functions 

^i(s), ro„(5), md Ot^Js) 

Because no concise expressions of the required func- 
tions are known, approximate formulas must be devised. 
The function 6^* has a simple operational equivalent, 
namely, 

D 



1 



(26) 



and, since the curves to be fitted are asymptotic in 
character, series of such functions were chosen as 
follows: 

ro«,(s)=5.75-3.756-^'2^Q^-1.50e~^-^^°* 

tf;r(s)^^3=0,096-0.053e-o-28i* 
?^r(s)4=6=0.045-0.032e~o.203s 

Figures 7 and 8 show the degree of exactness attained 
with these expressions. It was considered not im- 
portant to fit the curves accurately near the ongin. 



(27) 
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3, ha/ f chords 
Figure 8.— Growth of cir(njlatipn in two-dimensional flow, 

The operational equivalents To«,(-D), Wiw(P), etc., 
are easily written down from (26). The substitution 
of these equivalents into equations (22) and (25) and 
the evaluation of the resulting operators by the Heavi- 
side expansion theorem are quite lengthy and will not be 
reproduced. The resulting expressions for Ct^is) 
were foimd to be 



+ 0.043e-2-^^^^+0.915e-^'2'^^^ cos 0.095(5-19.135)] 
Ci^^^^ ir[l .589 - 0.242e-« ~ QAOSe~^ -^oo^ 
+0.008e-^ .m.+0.g72e-o.234s q 0734 (5-21 .1 17)1 



(28) 



Because the curves given by these formulas are 
considered in vahd near the start of the motion, new 
curves having the correct starting values given by 
equation (7) were drawn in as shown (fig. 9). These 
final curves have the useful approximate expressions: 



a£„^(s)^1.20e7r(l .000-0.283e-»-^0 
(7. (s)^L487r(1.000-0.36Ie-'>-38i*) 



(29) 
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An analogous expression for infinite aspect ratio is 
(7io^(s);:^27r(1.000-0.165e~o«*s*— 0,3^^^^^ (30) 
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s, half chords 
FiGUBE 9.— Indidal lift functions, ChJ^^ and CLais), 

LIFT IN VAKIABLE MOTION 

In addition to the lift given by the Hft function 
^Lai^)) airfoil experiences a reaction equal to the 
instantaneous rate of change of the normal-velocity 
component times the virtual additional mass of the 
wing in normal motion. In coefficient form: 



^ da 
Eds 



(31) 



Furthermore, if the wing is rotating in pitch, the effect 
of an additional relative camber is introduced. A 
simple integration, making use of weU-known results 
of thin-aixfoil theory, shows 

where the factor I is }i for a straight rectangular wing. 

1 TT 

For the elliptic wing, ^ ^^^ g' approximately, bemg 

somewhat smaller than % because the rotation intro- 
duces a smaller relative camber at the narrower sections 
toward the tip. 

The effects of combined vertical motion 

and rotation (a— ^) may be conveniently treated by the 
use of moving axes as shown in figure 10. With these 
points in mind, tlie following operational formula for 
the total lit fc may be derived: 

• OM =~Da{s) +d:,JD){ais) +lDB(s)] (33) 



w 

Figure iO.— Moving axes, a=w/C/ii, 



LIFT FUNCTIONS FOR AN OSCILLATING AIRFOIL 

The lift in sinusoidal motion where 



is given by 
Since 



... D 



(34) 
(35) 



1 



D 



(36) 



Expansion of this operator gives, with the exception of 
transient terms, 



(37) 



The function Cx„ {in) corresponds to the Hft func- 
tion G{n) introduced by Theodorsen (reference 2) for 
the oscillating two-dimensional airfoil, that is, in 
Theodorsen's terminology 

Oi^(in) =27rO(n) =2Tr{F{n) ~^iQ{n)] (38) 

The expressions for found from the operational 

eqtiivalents of (29) are: 



(^•+iGD..a=0.740-0.267^^:pO^ 



(39) 



Figure 11 shows these functions plotted against the 
wave ler^th 2t/7i of the oscillation. 



/in 
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Figure ll.—OseUialai^-lift functions, CLjin^^tviF^iG) and 'Ci.-(i«)=P+»Q. 

Kelation (37) is especially interesting (see reference 
7) because it shows a connection between the Fourier 
and the operational analyses. Thus, if the r^ponse of 
a linear system to a continuous sinusoidal excitation is 
known, 

Rn(s)=J{in)e'-^ (40) 
then, the function / immediately furnishes the opera- 
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tional equivalent of the unit response so that, for any 
yaiiable excitation s(s), 



Bis)^f(D)z(s)^JiDMD)l 



(41) 



LIFT IN GUSTS 



The foregoiag calculations provide the basis for the 
determination of lift under any prescribed conditions 
of motion of the wing. These results may also be used 
in conjunction with the equations given by Theodorsen 
(reference 2) to predict the air forces on wings with 
hinged flaps. 

In all cases treated, the airfoil has been considered 
as moving in air that would otherwise be at rest. An- 
other problem of considerable interest is the prediction 
of lift during passage of the airfoil through gusts. 
The two-dimensional theory for this case was developed 
by Kussner (reference 3) and has since been corrected in 
certain details by von Kdrmdn and Sears (reference 8). 

The basic solution required ia the gust problem is 
the solution for a unit sharp-edge gust of uniform up- 
ward velocity. In order to obtain this solution, it is 
useful to substitute for the change in direction of the 
relative air velocity an equivalent fictitious bending of 
the airfoil in still air such that it has at every point an 
angle of attack equal to the angle that would otherwise 
be produced by the gust. 

The effect of a bend progressing along the chord of 
the airfoil may be calculated by thin-airfoU theory 
(reference 9, chs. Ill and IV). A part of the effect 
appears as a change in angle of attack of the airfoil as 
a whole, namely: 



(42) 



The corresponding part of the hf t is obtained from the 
indicial-lift function Oi,^{s) by superposition. In 
addition, a reaction caused by acceleration of the non- 
circulatory potential flow exists during the time the 
airfoil is partly immersed in the gust. In two-dimen- 
sional flow, the additional reaction is 



AC.=2Vs(2-s) 



(43) 



No corresponding expression for the finite wing is 
known, but it may be reasoned that the maximum cor- 
rection will be no greater than that indicated by the 
inertia factor of the rigid elliptic disk, 1/E. Hence, the 
formula 



(44) 



was used for the finite wings as an approximation. 

The consideration of wings with curvature or sweep- 
back introduces another difficulty into the analysis, 
since the sections of such wings will not strike the edge 
of the gust simultaneously. It is obviously impractical 
to attempt to include in the analysis the effects of the 
many possible variations of plan form, arid the calcula- 

180817—40—2 



tions were therefore made on the assumption that all 
sections entered the gust simultaneously. Such an 
analysis may be considered sufficiently exact for the 
usual variations of plan form but is, of course, not ap- 
plicable to wings with considerable sweep. 

Figure 9 shows the functions, designated Cl^ (s), thus 
calculated. These curves have the useful approximate 
expressions: 

Ol,{^)a=z=1 .200?r(l .000"0.679e-*>-««8*--0.227e-3-2oO' 
<?i^(s)A=6 === 1.5007r(1.000---0.448e--<>*290s 

— 0.2726~»-'25*-0.193e-3-^00 (45) 
0/.^(s)A«c,=2.000ff(1.000-0.2366-o-^^«^ 

~ 0.5 1 36-^-2^^* 0, 1 7 1 e-2-^*> 
As in the case of the functions Ol^{s), the exponential 
forms were used to give simple operational equivalents. 
The operational equivalents of the indicial-gust func- 
tions, Gzffj give directly functions determining the al- 
temating lift of a stationary wing in an oscillating air 
stream. Thus 

CL(s) ='d^^ (m)^*«*= [P(n) +iQin)W^' (46) 

Figure 11 shows these functions in comparison with the 
corresponding fimctions for the oscillating airfoil. 




O 2 4 e a 10 12 i4 
s , half chords 

FxGTTBB 12.— Variation of the lift during passage through unit sharp-edge gust. . 



MOTION OF AIRPLANE IN GUST 



In most problems that arise in practice, the motion 
of the airfoil, or airplane, will not be prescribed before- 
hand but must be determined from dynamical equa- 
tions. The rising motion of an airplane (or, as it shall 
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be considered here, a loaded wing) while entering a 
sharp-edge gust presents such a problem and will be used 
to illustrate the application of the operational formulas. 

The dynamical equation for this case (neglecting 
pitching motion) is: 

m^+resisting force==iinpressed force (47) 

where the impressed force is that part of the lift caused 
by the gust. Since 

^^U± da 
dt '^cl2 ds 



duo 



2m p c Uq da_ £j'j 2^ 
gj^^2 27r2ds^^^2^rds 
2 



(48) 



2m 



where jLt=— • In coefficient form, 



(49) 



(5 
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F{0tTB3B 13.— M83dmiim-lift inorem^te developed In flying throogh a Unit 8liarp*edge 

gust. 



where ag is the change in angle of attack represented by 
the gust. 

For a unit sharp-edge gust, ag==l] then (solving 
for a), 

''""-Mm' 

By the use of the approximate expressions given for 
Cx^and (7i^ (equations (29) and (45)), this operator may 
be reduced to the form (19). 

Figure 12 shows the lift coefiBcient GLis)~(iDa{s) 
computed from equation (50) for several values of the 
density ratio /i and for J. ==6. Figure 13, derived from 
mmilar calculations, gives maximum lift loads attained 
in the sharp-edge gust as functions of the relative 
density. ' ' 
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iliiiiiipiiiiii^^^^^^ 

Po^ive direciipns of axes and angles (forces and naoments) are shown by arrows" 





(parall^ 
to axis) 
.gyinbol 


Moment about axis ^ ^ 


4^gie 


Yelpcities - 


Demgimtio^ 


iiiil 


Designation 


Bym- 

iiiif 


direction 


Designa- 


Sym- 

lliiii 


iliilKipiii 

(cdmpo-^ 
B€^taldng 
axis) 


AtJgular 


Lon^tudinal 

Lateral 


Y 
Z 


X \ 

Y 

Z 


BoUing 

PitcMng-^— 
Yawing 


L 

N 


Y >Z 

z — >x 

1 >Y 


Eon 

Pitch 

. Yaw 


0 


¥ 

V 

w 


P 

r 







Absolute coefliqienta Of moment 
(rolling) , (pitching) ; 



. . \ Angle of ^et of control surface (relative to neutral 

iV . positioji), d. (Indicate surface by proper subscript.) 

illllllllplll^ 



Diameter 

Geometric pij^cb \ ; . 
pJDf Pitch ratio 
F', Inflow velocity ^ . 

F^, Slipstream velocity \ \ 

3P, Thrust, absolute jC^o^ffioieat Cr= 



Torque, absolute epeflSicient €^=^ ^^5 



iBIiiiiliillllliB 

Q,^ S|)eed-power coefficientf^y ^ 

llllilllil 

lijplllillllH 



1 hp.=76.04 kg-m/s=550 ft-lb^/sec. 
i metric horsepower= 1 .013^ hp* 
i m.p.h. =0.4470 Da>p.3. 
1 m.p.s. =2.2369 m.pk. r ; 



^ 1 lb.=o!4536 kg. 

1 kg*^2,2046 lb. 
/ , ,1 jm.^l,609.35 m==6,2$0 ft. 
i>i=3,2808ft. 



